Locations of multicritical points for spin glasses on regular lattices 
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We present an analysis leading to precise locations of the multicritical points for spin glasses on 
regular lattices. The conventional technique for determination of the location of the multicritical 
point was previously derived using a hypothesis emerging from duality and the replica method. In 
the present study, we propose a systematic technique, by an improved technique, giving more precise 
locations of the multicritical points on the square, triangular, and hexagonal lattices by carefully 
examining relationship between two partition functions related with each other by the duality. We 
can find that the multicritical points of the ±J Ising model are located at pc = 0.890813 on the 
square lattice, where Pc means the probability of Jij = J(> 0), at Pc ~ 0.835985 on the triangular 
lattice, and at pc = 0.932593 on the hexagonal lattice. These results are in excellent agreement with 
recent numerical estimations. 



I. INTRODUCTION 

The realistic world is affected by randomness. Unlike 
non-random systems, those with randomness sometimes 
show rich and complicated phenomena. One of the inter- 
esting issues for such random systems is spin glass. 

Many studies mainly by the mean-field analysis have 
been successful to elucidate various concepts for under- 
standing spin glasses [H, 0, Q • One of the current issues 
in spin glasses is their nature in finite dimensions below 
the upper critical dimension. Unfortunately, for finite 
dimensions, we often rely on numerical simulations, be- 
cause there are few ways to analytically study spin glasses 
in finite dimensions. We need long equilibration times 
for the numerical simulations for spin glasses and aver- 
age over many realizations of random systems to make 
error bars small enough. It is thus difficult to give con- 
clusive understanding on nature of spin glasses in finite 
dimensions. 

To establish reliable analytical theories of spin glasses 
has been one of the most challenging problems for years. 
A part of successful analyses to elucidate properties on 
spin glasses is by the use of the gauge symmetry. By use 
of the gauge symmetry, one can obtain the exact value 
of the internal energy, evaluate the upper bound for the 
specific heat, and obtain some correlation inequalities in 
a subspace known as the Nishimori line [1, Q . This gauge 
symmetry also enables us to rewrite the free energy along 
the Nishimori line as the entropy for the distribution of 
frustration Q. Many aspects on spin glasses are essen- 
tially related with frustration. Therefore we expect the 
possibility that a basis of establishment of a systematic 
approach to spin glasses would be in the gauge symmetry. 

A recent related development with the gauge symmetry 
is the conjecture to predict the location of the multicrit- 
ical point, which is the special point lying on the inter- 
section between phase boundaries and the Nishimori line 
as in Fig. [J] 0, i, i, [HI, d • The predictions by the 
conjecture have shown agreement with numerical estima- 
tions roughly with precision to the third digit. The con- 
jecture has opened a way of a general scheme for determi- 



FIG. 1: Phase diagram of the ±J Ising model on two- 
dimensional lattice (left panel) and on higher dimensions 
(right panel). The vertical axis expresses the temperature 
T, and the horizontal line denotes the concentration p of the 
antiferromagnetic interactions. The multicritical point is de- 
scribed by the black point (MCP). The Nishimori line is de- 
scribed by the dashed line. For higher dimensions, not only 
the ferromagnetic (FM) and paramagnetic phases (PM) but 
also the spin glass phase (SG) exists. 



nation of the location of the multicritical points for spin 
glasses on any self-dual lattices and mutually dual pairs 
of lattices. Nevertheless it has been found that the con- 
jecture on several mutually dual pairs of hierarchical lat- 
tices does not always give predictions in agreement with 
estimations by the renormalization group analysis p^ . 
Such discrepancies are not negligible, because the renor- 
malization group analysis on hierarchical lattices gives 
exact solutions. To construct a more reliable technique, 
we have improved the technique leading to the location 
of the multicritical point by combining the concepts of 
the renormalization group analysis with the duality [3] . 
The improvement for hierarchical lattices has greatly suc- 
ceeded as seen in the literature, because the discrepan- 
cies between the predictions by the improved technique 
and the exact estimations by the renormalization group 
analysis actually decrease. Reconsidering the improved 



2 



technique on the hierarchical lattices, in this paper, we 
apply the improved technique to the regular lattices such 
as the square, triangular, and hexagonal lattices. 

The present paper is organized as follows. In Sec. II, 
we introduce the conventional conjecture and take a look 
at several published predictions. In addition, the prob- 
lem on the conjecture is pointed out here. The improved 
technique on the regular lattices is proposed after review 
of the case on the hierarchical lattices in Sec. Ill, and 
formulated in Sec. IV. The improved version shows a 
very close relationship with the entropy of the distribu- 
tion of frustration as shown in this section. In Sec. V, we 
carry out the explicit calculations by the improved tech- 
nique for the regular lattices. Moreover we have to care- 
fully evaluate the performance of the improved technique, 
comparing their predictions with the existing results. We 
examine the correspondence with the Domany's exact re- 
sult of the slope of the critical point on the phase 
diagram. The conclusion is in the last section of the 
present paper. Sec. VII. 



II. CONVENTIONAL CONJECTURE 
A. Duality 

It will be useful to review the analysis by the duality, 
for the conjecture is established by combination of the 
duality and the replica method 7, 8]. The duality is one 
of the tools to identify the transition points for various 
types of classical spin systems such as the Ising model, 
and the Potts model by use of a symmetry embedded 
in the partition function [13]. We take the non-random 
Ising model on the square lattice as an example. The 
partition function is given as 



{S.} fe> 



(1) 



Where Si is the Ising spin taking ±1 and the prod- 
uct with the subscript (ij) is over the nearest neighbor- 
ing sites. We can regard this partition function as the 
multi-variable function of components of the edge Boltz- 
mann factor. In this case, the edge Boltzmann factor is 
exp {f3JSiSj)). We consider that the dependence on (3 
emerges through this edge Boltzmann factors. We set 
two components as xo{(3) = e^'^ and xi{j3) = e~^'^ for 
convenience. The component xo{(3) is often called the 
principal Boltzmann factor, which is defined by the edge 
Boltzmann factor for the state with all edge spins parallel 
The principal Boltzmann factor is an important compo- 
nent throughout this paper. The duality is carried out 
by the Fourier transformation for this edge Boltzmann 
factor defined on each bond of the lattice [l^. The two- 
component Fourier transformation gives the dual edge 



Boltzmann factors as. 
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,/3J , 



(2) 
(3) 



As a result, we establish the relation between the parti- 
tion functions with different components as. 



Z{xo,xi) = Z{xq,xI). 



(4) 



We here extract two principal Boltzmann factors xq and 
Xq to measure the energy from the state with edge spins 
on each bond being parallel as. 



x^-{(3)z{u,)^xl"-{f3)z{ul) 



(5) 



where Nb stands for the number of bonds, and ui and 
ul are called the relative Boltzmann factors defined as 
= xi{P)/xo{P) and ul{(3) = xl{(3)/x^{P). Each 
partition function is now reduced to a single-variable 
hmction of ui and u*, whose explicit forms are. 



MP) 



-2/3J 



tanh pj. 



(6) 
(7) 



Being very well known, the duality relation can be given 
as e~'^^ = tanh/? J to think of the dual partition func- 
tion as one of another Ising model with the edge Boltz- 
mann factor qP'J'^^'', , We can identify the critical point 
as a fixed point of the duality as e^^''^ = tanh /3c J, under 
the assumption of a unique transition. On this critical 
point, an appealing equation is satisfied xq = Xq. 



B. Duality for the quenched system 

We consider the case of the random-bond Ising model 
and review the conventional conjecture . The Hamil- 
tonian is defined by 



{ij 



ijSiSj, 



(8) 



where denotes the quenched random coupling and the 
summation is over the nearest neighboring sites. Though 
various types of distribution for can be considered, we 
here restrict ourselves to the ± J Ising model for conve- 
nience. The distribution function for the ± J Ising model 
is given by 



P(Jy) = p6iJ,j-J) + il-p)SiJ,j+J) 



2 cosh /3p J 



, (9) 



where /3p is defined by e~'^^p'^ = {l—p)/p. The Nishimori 
line is given by the condition (3 = (3p and is described by 
the dashed line in each phase diagram of Fig. [T] 



We apply the replica method to the ±J Ising model 
on the Nishimori line as, 



^(/3p,/3p) = 



J2 n ^MPpJvsrsf) 

{S,}"=1 



(10) 



where n stands for the replica number and the angu- 
lar brackets denote the configurational average. We ap- 
ply the duality argument as reviewed above to this n- 
replicated ±J Ising model. The duality gives the rela- 
tionship of the partition functions with different compo- 
nents of the edge Boltzmann factor, 



Zn {Xq , Xi, ■ ■ ■ , X„), 



(11) 



where the subscript of x and x* denotes the number of 
antiparallel-spin pairs among the n replicas. Two princi- 
pal Boltzmann factors xq and Xq are given as 0, [11 , 
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(12) 
(13) 



We extract these principal Boltzmann factors similarly 
to the case of the non-random bond Ising model. 



Xo{f3)^'' Zn{ui,U2, ■ ■■ ,Un) 



(14) 



We remark that the partition function for the n- 
replicated ±J Ising model is a multi- variable function 
of the edge Boltzmann factors yet differently from the 
case of the non-random bond Ising model. 



C. Conjecture 

We describe schematically the relationship two re- 
duced partition functions the curves of the rel- 
ative Boltzmann factors (ui(/3), U2(/3), • • • ,Un(/3)) (the 
thin curve going through the multicritical point pc) and 
(wJ(/3),M^(/3),- • • ,<(/?)) (the dashed line) as in Fig. H 
We now consider the relationship between these curves by 
the projections on the two-dimensional plane (ui, U2) for 
convenience. As the temperature changes from to cx3, 
the representative point (ui(/3), M2(/3), • • • ,«„(/?)) moves 
toward the point P (the high-temperature limit) along 
the thin line in Fig. [21 Then the corresponding dual 
point (uj(/?), U2(/3), ■ • ■ ,u*(/3)) moves along the dashed 
line in the opposite direction from P to F (the low- 
temperature limit). These features have been shown rig- 
orously and imply the existence of the duality relation for 
the temperature [l2| • If two curves describing change of 
(mi(/3),U2(/3),-- - ,M„(/3))_and (ut(/?), (/?),•• • ,<(/3)) 
become completely coincident with each other, we can 
obtain a relation Ur{l3*) = u*((3). Solving this relation, 
we obtain the duality relation for the temperature (3* {(3) . 




FIG. 2: Schematic picture of the duality and the renormal- 
ization flow for the replicated ±J Ising model. 



The well-known duality relation e~^^ — tanh/3J for 
the non-random Ising model can be indeed derived from 
the relation Ur{P*) = u*{(3). Unfortunately the thin 
curve (iti (/?) , M2 (/3) , • • • ,u„(/?)) does not coincide with 
the dashed curve {u\{l3), u^iP), • ■ ■ , u^{l3)) for the repli- 
cated ±J Ising model on the Nishimori line as shown 
schematically in Fig. [51 In the case for the replicated 
±J Ising model, we can neither find a duality relation 
for the temperature explicitly, nor identify the multicrit- 
ical point as the fixed point of duality. 

We thus have provided a hypothesis on determination 
of the multicritical point 0, [1]. We assume that the 
equation xo(/3) = Xq(/3) is also satisfied at the multicriti- 
cal point, similarly to the non-random bond Ising model 
at the critical point, for any n-replicated systems includ- 
ing the quenched system (n ^0) 0, H, i, 0, HH, [ll- 
Validity of this hypothesis can be rigorously shown for 
n = 1 and 2, and for n = 00, and has been numerically 
confirmed for the n = 3 case of the replicated ± J Ising 
model on the square lattice within error bars of the nu- 
merical simulation % Relying on these facts, we assume 
that we can predict the location of the multicritical point 
for the ± J Ising model on the square lattice by the single 
equation xq[I3) — Xq{B). The quenched limit rt — > for 

3 0,i, 



this equation yields 



1 



plogp- (1 -p)log(l -p) = 2 log 2- 



(15) 



The solution to this equation is pc — 0.889972. We can 
also predict the location of the multicritical point for 
other types of randomness. For instance, that of the 
Gaussian Ising model with the mean Jo and the variance 
J = 1 is given as Jq = 1.021770 Si]. 

The conjecture can be extended to mutually dual pair 
lattices by considering the product of two partition func- 
tions with different couplings indicating two multicrit- 
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type conjecture numerical result 



SQ ±J 


Pc 


= 0.889972[7, 8] 


0.8905(5) [19] 








0.8906(2) [20, 21] 








0.8907(2) [22] 








0.8894(9) [23] 








0.8900(5) [24] 








0.89081(7) [25] 


SQ Gaussian 


Jo 


= 1.021770 [7,8] 1.02098(4)[21] 


TR ±J 


Pc 


= 0.835806[11] 


0.8355(5) [24] 


TR Gaussian 


Jo 


= 0.798174 




HEX ±J 




= 0.932704[11] 


0.9325(5) [24] 


HEX Gaussian Jo 


= 1.270615 





TABLE I: Comparisons among the results derived by the con- 
ventional conjecture and the existing numerical results. SQ 
denotes the square lattice, TR means the triangular lattice, 
and HEX means the hexagonal lattice. 



ical points jlOj . and we can obtain a relation between 
these multicritical points. We can predict the location of 
the multicritical point on the triangular lattice, though 
it is the mutually dual pair with the hexagonal lattice, 
by simultaneous use of the star-triangle transformation 
Through these extensions of the conjecture, we have 
predicted locations of the multicritical points for various 
cases, which show good consistencies with other estima- 
tions with high precision in spite of the simplicity, as seen 
in Table U 

On the other hand, the conjecture has given approxi- 
mate locations of the multicritical points on several hier- 
archical lattices with the largest discrepancy of 3% [l^ . 
We have inferred that the issue yielding such discrep- 
ancies comes from the following fact. Unlike the non- 
random Ising model, two lines do not coincide as shown in 
Fig. [2] Nevertheless we have assumed that the equation 
xo = Xq is satisfied. As shown in next section, the tech- 
nique used in the conjecture has indeed been improved 
by reconsideration of these problems. 



III. IMPROVED TECHNIQUE 

A. Renormalization Flow 

We here introduce a new point of view from renor- 
malization group, which has helped us to improve the 
technique on the location of the multicritical point on 
the hierarchical lattices In the present study, we 

apply this idea to the regular lattices. Two features 
through the renormalization group analysis on the hi- 
erarchical lattices are found, (i) The partition function 
does not change its functional form by renormalization; 
only the values of coupling constants change because of 
specialty of hierarchical lattices, (ii) The renormalization 
flow starting from a critical point is attracted toward an 
unstable fixed point. The feature (i) permits us to ex- 



press the renormalization flow following the arrows ema- 
nating from Pc and dc to C, ph and di to P, and pi and dh 
to F as in Fig. O By the feature (ii), there is the renor- 
malization flow from the multicritical point pc and corre- 
sponding dual point dc, which reaches the unstable fixed 
point C, {u^\u^2°\--- where the superscripts 

mean the number of steps of renormalization. Consid- 
ering the two properties of the renormalization flow, we 
can find that the duality relates two trajectories from pc 
and from dc, tracing the renormalization flows at each 
renormalization. In other words, after a sufficient num- 
ber of renormalization steps, the thin curve representing 
the original system and the dashed curve for the dual sys- 
tem both approach a common renormalized system rep- 
resented by the bold line in Fig. [21 which goes through 
the unstable fixed point C. On this bold line, the par- 
tition function can behave as a single-variable function 
and we can identify the unstable fixed point by a sin- 
gle equation xo{K) = Xq{K) similarly to the case of the 
non-random Ising model. This fact enables us to assume 
to predict the exact location of the multicritical point by 
the following equation, 

xt\P)=4^°^\P). (16) 

However it is difficult in general to evaluate this equa- 
tion. We therefore have proposed a first-approximation 
equation Xq^^ = Xq*-^^ as the improved technique. This 
equation has indeed led to more precise results than the 
relation xo{(3) = Xq(/3) does |l4j |. 

B. Partial Summation 

We step in the stage of establishment of the improved 
technique on regular lattices. Unfortunately, if we apply 
the above renormalization group analysis to the regular 
lattices, the feature (i) may be incorrect, because other 
types of many-body interactions are generated after each 
renormalization. In addition, the renormalization group 
analysis for the regular lattice is usually regarded as an 
approximate tool, because many-body interactions ap- 
pear after the renormalization which prevents us from 
iterating the renormalization. If we attempt to construct 
the recursion relation of the renormalization for regular 
lattices, we introduce some approximations. However we 
recall that the improved technique on the hierarchical lat- 
tice has been successful within satisfactory precision even 
by a one-step renormalization 14]. We then sum over in- 
ternal sites only in a unit cell of each hierarchical lattice. 
For example, we show such a calculation explicitly for 
the non-random Ising model on a self-dual hierarchical 
lattice in Fig. [31 

Sl,S2 

(17) 

where A is the extra coefficient yielded by the renormal- 
ization. The left-hand side corresponds to the principal 
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r=0 r=l r=2 

/3(2) ^(1) ^(0) 



FIG. 3: A self-dual hierarchical lattice. The number r repre- 
sents the step of the construction of the hierarchical lattice. 
The renormalized inverse temperature is expressed by /?''', 
where I is the step of the renormalization. 

Boltzmann factor after one-step renormalization when 
SiSj = 1. We can regard this renormaUzed principal 
Boltzmann factor as the partition function defined on 
the unit cell of the hierarchical lattice in Fig. [3] under 
the constraint that the spins at both ends are parallel. 
In other words, the principal Boltzmann factor for the 
improved technique is in general given as 

unit 

where the overline means the summation over internal 
spins and ^2 on the unit cell of the hierarchical lattice 
with the spins on the edge up Si and Sj = 1. The prod- 
uct runs over the nearest neighboring pairs on the unit 
of the hierarchical lattice. For the case for the random- 
bond Ising model, after application of the replica method, 
we also obtain the similar principal Boltzmann factor to 
the above example, by performing the partial summation 
over the unit of the hierarchical lattice. Instead of iterat- 
ing summation as the renormalization, we sum partially 
over sites on the regular lattice. After we partially trace 
out the degrees of freedom in some area of the regular lat- 
tice, we regard the partition function on this area as the 
principal Boltzmann factor as the above example. We are 
able to construct the improved technique for the regular 
lattice without neglecting many-body interactions emerg- 
ing after the summation. We call the limited-summation 
area of the regular lattice the cluster in this paper, and 
we define the partition function on the cluster as the prin- 
cipal Boltzmann factor for the improved technique on the 
regular lattice. When we take the summation over the 
internal sites, we impose the fixed boundary condition 
on all the spins on the boundary of the cluster as we 
have done for edge Boltzmann factor in the conventional 
conjecture. 

We remark the possibility of the performance of the im- 
proved technique. Even if we partially sum over degrees 



of freedom on the regular lattice as considered above, it 
is necessary to use an infinite-dimensional space to si- 
multaneously express changes and generations of various 
types of many-body interactions. In this infinite dimen- 
sional space, we assume that there are two renormal- 
ization flows going uniformly toward the unstable fixed 
point, similarly to the case on hierarchical lattices. It is 
not able to examine such behaviors of the renormaliza- 
tion flows in the infinite-dimensional space and we cannot 
verify this assumption. Equation xq = Xq on the cluster 
may hence be a worse approximation than the conven- 
tional conjecture x'^^ = Xq^^^ by consideration on a single 
bond. We have to carefully evaluate the performance of 
the improved technique for the regular lattices. At least, 
the obtained results as seen later can, however, give an 
answer for the location of the multicritical point for the 
± J Ising model on the square lattice with the precision 
to the fourth digit. 

In addition, we will consider several types of the im- 
provements for the regular lattice. On the hierarchical 
lattices, sufficient renormalization enables us to obtain 
the exact solution. Therefore we can expect that the im- 
proved technique for the hierarchical lattices gives more 
precise answers if we use the principal Boltzmann factor 
after more renormalization steps. On the regular lattices, 
it is considered that the performance of the improved 
technique depends on the area of the cluster. The num- 
ber of degrees of freedom partially summed over on the 
regular lattice would correspond the step of the renor- 
malization on the hierarchical lattice. We propose a few 
types of the clusters including several bonds below, con- 
sidering this assumption. 



IV. FORMULATION 

A. Square Lattice 

The starting point for the establishment of the im- 
proved technique is the exact duality relation (jlip for 
the n-replicated partition function. As shown in Fig. 31 
we consider to sum over a part of the spins on the square 
lattice. Then the exact duality relation (fTTj) is reduced 
to 

'7(s)('^(«) ™(«) ... r(^h - Zt^^rr*'"^ r*^^) ••. t*(*)1 
\-^0 1-^1 1 T^n ) ~ V-^O '•''1 ' T-'^n )■ 

(19) 

Here Zn^ represents the reduced partition function by 
the summation of a part of spins on the square lattice. 
The superscript s distinguishes the type of the approx- 
imations, for we consider different types of the summa- 
tion below. The quantity x\ is the edge Boltzmann 
factor including many-body interactions generated after 
the summation. We take a cluster of the square lattice 
as in Fig. [?] and define the principal Boltzmann factors 
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FIG. 4: Example of the partial summation for the square 
lattice and the cluster. The top figures describe one of the 
types of the summation for the square lattice. The bottom 
figures express the cluster for the principal Boltzmann factor. 



after the summation, 



part 

part 



(20) 



-0Ji 



V2 



SiSj) 



(21) 



where the overline means the summation over internal 
spins in the cluster of the square lattice as the filled cir- 
cle in Fig. IHwith the other spins fixed to up directions 
{Si = 1} represented by the white circles in Fig. [4l The 
word "part" represents that the product runs over the 
bonds of the cluster under consideration. These princi- 
pal Boltzmann factors can be regarded as partition func- 
tions after the configurational average and application of 
the replica method defined on the cluster under the fixed 
boundary condition as shown in Fig. 3] 

We assume that a single equation gives the critical 
points for any number of n, similarly to the conventional 
conjecture, 



(22) 



By the extrapolation of the quenched limit n ^ of this 
equation, we establish the improved technique for the 
square lattice as follows, 



logZ*(^)(/3,{J,,})l - [logZ(^)(/3,{J,,}) 



= 0. 

(23) 



We need the configurational average for of the loga- 
rithmic terms by two partition functions Z*^"-* and Z**^*^ 



defined on the cluster, 



part 



Z(^)(/?,{J,,}) = ^ J^e^''-^^-^^- (24) 

part 

{S.} (ij) ^ 

(25) 

where the asterisk means that the edge Boltzmann factor 
is given in a different form obtained after the duality for 
the ± J Ising model. We can estimate the location of the 
multicritical point by the above relation ([23)) as detailed 
below. 



B. Triangular Lattice 

Similarly to the case on the square lattice, we can de- 
rive the improved technique for the triangular lattice. 
We give first several remarks on the duality with the 
star-triangle transformation. For the triangular lattice, 
it is convenient to use the face Boltzmann factor instead 
of the edge Boltzmann factor. The original face Boltz- 
mann factor has three edge Boltzmann factors defined 
on the bonds of the elementary triangle. On the other 
hand, the dual face Boltzmann factor has three dual edge 
Boltzmann factors defined on the bonds on the star shape 
as in Fig. [5] [HI, H^. In the dual face Boltzmann factor, 
the summation over the spin at the center of the star is 
included, which corresponds to the star-triangle trans- 
formation. Similarly to the case of the square lattice, we 
start from the exact duality relation for the replicated 
random-bond Ising model on the triangular lattice. 



Z' 



TR,n 



(^0,^1,- 



Z^'^l (A* A* 



(26) 



where A and A* are the original and dual face Boltzmann 
factors. We consider the summation over a part of the 
spins on the triangular lattice as in Fig. [51 and then this 
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FIG. 6: Example of the summation for the triangular lattice 
and the cluster. The top figures express one of the types of 
the summation for the triangular lattice. The bottom figures 
express the cluster for the evaluation of the principal Boltz- 
mann factor. 



equation is reduced to, 



(Ay") ^i*' 



■ ) _ z^1{Aq 



(s) 



(27) 



(s) 

where Z!^^ ^ represents the reduced partition function 
by the summation over a part of spins on the triangular 

(s) 

lattice. The quantity A), is the face Boltzmann factor 
after the summation. We define the principal Boltzmann 
factor after the summation over a part of the triangular 
lattices as, 



A, 



*{s) 



part 



g/3Jl2SlS2+/3J23S2S3+/3J3l5'3Sl 

{S,} A 



(28) 



part 



EEn 

{S.}{So} A 



V2 



Si So) 



1 



V2 



fi.ho I „-/3J20 



S2Sq) 



X [e'^--^" + e 



(29) 



where ivj'^ is equal to the number of the up-pointing 
triangles included in the cluster. The overline means the 
summation over internal spins included in the cluster of 
the triangular lattice as the filled circles in Fig. [Hlwith 
the other spins (the white circles) up {Si = 1}. The 
word "part" represents summation over the up-pointing 



triangles in the cluster as in Fig. [6l The summation 
over 5*0 means the star-triangle transformation. We can 
establish the improved technique for the triangular lattice 
by an equation of these principal Boltzmann factors after 
the summation A\^^ — A*^^^^^ and by the extrapolation to 
the quenched limit of n — > 0, 



logZ;(^)(/3,{J,,}) 



logZ(f^(/3,{J,,}) 



(30) 

where the partition functions Z'!^^ and .^■^'j^'' on the clus- 
ter of the triangular lattice are defined as. 



_ part 



g/3Jl2SlS2+/3./23S2S3+/3J3lS3Sl 

{S.} A 



(31) 



part 



EEn 

{S.}{So} A 



V2 



1 I 

X —= (e 
V2 ^ 
1 



PJ^iO _j_ g — /3J20 



5250) 



7i<' 



S3S0) 



e-^''°SiSo) 



(32) 



C. Frustration Entropy and Multicritical Point 

Before going into the detailed calculations for deter- 
mination of the locations of the multicritical points, we 
consider the physical meaning of the improved technique 
from a different point of view. We show the relationship 
between the improved technique and the gauge symme- 
try, by taking the case of the square lattice as an ex- 
ample. The observation below implies the existence of 
deeply physical meaning behind the multicritical point. 
We will discuss the structure of the phase diagram by 
use of the connection between the improved technique 
and the gauge symmetry as shown below in the following 
section. 

The second term Z^'^^ on the left-hand side of Eq. ([23)1 
can be regarded as the free energy for the random-bond 
Ising model defined on the cluster if divided by /?. The 
free energy for the gauge invariant model can be reduced 
to another form by the gauge transformation defined as 



J, 



S^ 



J ij ^ j 



(33) 
(34) 



where <Ji takes either —1 or -t-1. We can obtain a useful 
expression of log Z^'') for the case of the ± J Ising model 
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o 



a- 



6 



of Eq. ^ as 



\0gZ<^\(3,{J^J}) 



N 



f -1 log2. 



(38) 



FIG. 7; Dual lattice for the cluster in Fig. H The black 
colored sites are free spins as the targets of the summation 
and the white ones are fixed in the up-pointing directions. 



for instance 



logZ(^)(/3,{J.,}) 



The first term on the right-hand side of this relation can 
be reduced to the same expression as Eq. ([35)) . 



log {J,,}) 
1 



2^D'(2cosh/3pJ)^B' 

X J2 4'^(/3p,{J.,})log4''(A{>/..})- (39) 



2^vi"' r.«ot, « 7^JV<=' 



{2 cosh Ppjy 
X ^ z(^)(/3p,{J,,})logZ(^'(/3,{J.,}), (35) 



(s) 

where is the number of bonds included in the cluster. 
In this expression, we obtain the entropy of the distribu- 
tion of frustration on the cluster by setting f3 = (3p. 

On the other hand, the first term Z*^^'^ on the left- 
hand side of Eq. which is generated from the 
dual principal Boltzmann factor after the summation, 
is not in a gauge invariant form. The duality, how- 
ever, can transform Z**^*' into a gauge invariant form. 
At this time, the form of the edge Boltzmann fac- 
tor changes from {exp(/3Jy) + e-xjp{—(3Jij)SiSj} /\/2 to 
exp(/3Jy5,S'j) m. 



(36) 



where Ng'''^ is the number of sites in the cluster of the 
original square lattice. However the form of the lattice 
is transformed as shown in Fig. [71 We denote this fact 
by the subscript D. The apphcation of the duality to the 
partition function Z*^'') on the cluster with the dual edge 
Boltzmann factor enables us to derive a gauge invariant 
form, 



part(D) 

E n 

{S^} 



(37) 



where the 'part(D)' expresses the product over the bonds 
on the dual cluster as in Fig. [T] The configurational- 
averaged quantity of its logarithm can be regarded as the 
free energy of the random-bond Ising model defined on 
the dual lattice for the small lattice. We use the relation 
([55)1 and rewrite the first term Z*^*^ on the left-hand side 



(s) 

Here expresses the number of sites of the dual clus- 
ter, which is equal to that of the plaquettes of the original 
cluster. In the case of the cluster of the square lattice as 



in Fig. m ^i'^ = 1, 



4, N 



D 



considerations enable us to rewrite Eq. 



4. The above 
as. 



Nf^ + 1 log 2, 



(40) 



where 



{j.^} (2cosh/3pJ)^«' (2cosh/JJ) 



N 



(41) 



{j,^} (2 cosh /3p J) 



(2 cosh 13 J) 



Nl 



(42) 



If we set P = Pp in this expression, we find that Eq. ((40)) 
states that the multicritical point is located where the 
difference between two entropies of the distribution of 
frustration on two lattices related by the duality takes a 
special value. We can use this expression of the improved 
technique to lead the structure of the phase diagram such 
that the gauge symmetry and several correlation inequal- 
ities used to give the prediction as in Fig. m\^M- 

It is 

straightforward to establish the expression as in Eq. (jlQ]) 
also for the case of the triangular lattice. 

In the next section, we show the results for the precise 
locations of the multicritical points obtained by comput- 
ing Eq. for the square lattice and Eq. (150)) for the 
triangular lattice. 
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FIG. 8: Two patterns of the clusters for the improved tech- 
nique on the square lattice. The top figures express the clus- 
ters for the evaluations of Z^^^ and Z*'"'. The bottom figures 
represent the dual lattices for the clusters, on which the par- 
tition functions are denoted by The filled circles are 
the targets of the summation and white ones are fixed to up 
directions {Si} — 1. 



V. DERIVATIONS OF MULTICRITICAL 
POINTS 



We derive the location of the multicritical point by the 
improved technique for the regular lattices. If we con- 
sider a larger range of the summation of spins, (i. e., the 
cluster includes more bonds and sites.) it is expected that 
the precision of the improved technique becomes higher. 
One of the reasons is that the improved technique can in- 
clude more effects of spatially non-uniform interactions, 
which are essential features in random spin systems. In 
other words, the conventional conjecture has been the 
zeroth approximation without consideration of a form of 
the lattice and non-uniform interactions in space. The 
improved technique is also an approximation but gives 
more precise answers than the conventional conjecture, 
because it is formulated with the consideration of an in- 
dividual characteristic of the lattice similarly to the case 
of the hierarchical lattices. We express the type of the 
approximations by the value of s, which has represented 
the form of the cluster. In this paper, we show the deriva- 
tions of the multicritical point by use of several types of 
the clusters for the cases on the square lattice as in Fig. 
[5] and on the triangular lattice as in Fig. [51 




FIG. 9: Two patterns of the clusters for the improved tech- 
nique on the triangular lattice. The left-hand side for each 
type of the approximations is for the partition function .^^ri 
and the right-hand side is for the dual partition function ^^a' ' 
for which the star-triangle transformation is needed. We use 
the same symbols as in Fig. [S] 



A. First Approximation for the Square Lattice 

We consider the first approximation for the location of 
the multicritical point on the square lattice by the s = 1 
cluster as shown in Fig. [H To identify the multicritical 
point, we solve Eq. We first calculate the partition 

function on the cluster in Eq. (j24p . The summation over 
the single spin Sq at the center surrounded by four bonds 
(Jqi, Jq2, Jq3, Jqa) with four spins up yields 



E 



So=±l 

= 2 cosh{/3( Joi + Jo2 + Jo3 + ^04)}. (43) 
Another partition function in Eq. ([M)) is calculated as, 

So=±li=l 

4 r 4 



1 \ 

-= m(2cosh/3Jo.) + n(2sinh/3Jo. 
v2/ [j^i j=i 



(44) 



This quantity is also obtained from the evaluation of the 
partition function Z^"^ defined on the dual cluster by 
use of the relation ([M]) . We can explicitly write down 
the equation for the precise location of the multicritical 
point for the ±J Ising model on the square lattice as, 
fromEq. 

+ tanh^ Kp J[ To}j log + tanh^ K To}j 

I 2 cosh {KpJ^Loroi} 2 cosh {kJ^Lo ^o. } 
- o , ,4 log ■ ' 



2 (2coshis:p) 
2 log 2, 



(2 cosh X)' 



(45) 
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where we use the couphng constant K = (3 J and its sign 
Tij. Setting K = Kp, we solve this equation and ob- 
tain pi^'' = 0.890725. This result is listed in Table |TT] to 
see the performance of the improvement and to compare 
the results by the improved technique with the existing 
ones. Another type of the approximations for the square 
lattice is by the cluster labeled by s = 2 in Fig. [5] and 



J 



can be straightforwardly evaluated. The numerical ma- 
nipulation of this approximation give another prediction 

-- 0.890822 as listed in Table El 



For the Gaussian Ising model, we have to evaluate the 
quadruple integration over four bonds {Jij} for the im- 
proved technique as, 



n P{Jor)dJor log I n (2 cosh/3Jo,) + Y[{2 sinh/3Jo.) [ 

° i=l U=l 2=1 J 

/oo 4 
Yl P{.hi)djQ, log {2 COSh{/3( Joi + Jo2 + ^03 + Jo4)} = 2 log 2, 
'OO ■_i 



(46) 



where P{Jij) is the Gaussian distribution function with 
the mean Jq and the variance J — 1. The numerical 
manipulation of this equation gives the location of the 
multicritical point for the Gaussian Ising model on the 
square lattice as J^^^ = 1.021564. 



B. Other Approximations for the Square Lattice 

We restrict ourselves to the case of the square lattice 
and consider other types of the improvement. The key 
of the improvement is the cluster reflecting the shape 
of the square lattice. We consider here another type of 
approximations by dividing the square lattice into two 
clusters, which can cover the whole of the lattice, as in 
Fig. [ini Then we can reduce the duality relation pT|) to, 



1^) 



=1 : 



.^(t) At) 



= Zl:^^{x^^,^^,---,xf>,x^,---l (47) 



is t) • 

where Zn ' is the reduced partition function after the 
summation over internal sites on two clusters, and x^*-* 
and a;*^*-* are the edge Boltzmann factor including many- 
body interactions after the summation denoted by s and 
t, respectively. The reduced partition function is re- 
garded as a multi-variable function of two types of ar- 
guments a;'-"' and x^*\ We extract the principal Boltz- 
mann factors and assume that a single equation gives the 
location of the multicritical points. 




s = 1 



















) 











































/ 

















FIG. 10: Other approximations by use of two clusters on the 
square lattice. The small lattice t = \ includes 12 bonds and 
t = 2 has 24 bonds. 



{s) (t) 



*{s) *(t) 



(48) 



From this equation, we estimate the location of the mul- 
ticritical point. For s = 1 and t = 1 small lattices, we 



obtain p, 



(1,1) 



lattices, pi 



(1,2) 



0.890794 and, for s = 1 and < = 2 small 
= 0.890813. If the cluster includes many 



bonds, we expect that the improved technique can ap- 
proach the answer. It is then considered that the mul- 
ticritical point is located at = 0.890813. At least, we 
find that both of the estimations indicate a higher value 
Pc ~ 0.8908 than ■p'c^ k. 0.8900, and the precise location 
of the multicritical point would be pc ~ 0.8908. 
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C. First Approximation for the Potts Spin Glass 

As another application of the improved technique, let 
us consider the Potts spin glass on the square lattice de- 
fined by the Hamiltonian, 



H 



(49) 



where = 



- 4>j expresses the difference between ad- 
jacent Potts spins taking an integer value between and 
q — 1. The quantity lij is the random variable following 
the distribution function given as 



l-{q- l)p {k, = 0) 
P{Uo) = { P {hj ^ 0) 



e^p + q — 1 



,(50) 



We here estimate the location of the multicritical point 
for the Potts spin glass on the square lattice by use of the 
s = 1 small lattice. The partition functions on the cluster 
as in Eq. (^5]) are given as. 



(56) 



part 



part 



(57) 

, / (I \. 7! J 

We can carry out the summation over 00 at the center of 
the s = 1 cluster as in Fig. [8] as. 



where e^^ = {1 — (q — l)p}/p. This Potts spin glass also 
has the gauge symmetry. For the Potts spin variables and 
random variables, we define the gauge transformation as, 



0i + Si, 

hj ~ {Si — Sj). 



(51) 
(52) 



Here Si takes an integer between and q — I. Therefore 
we can establish the Nishimori line by setting f3J = Kp, 
where the internal energy can be calculated exactly and 
the specific heat can be bounded [l^]. The edge Boltz- 
mann factor for the Potts spin glass is given as. 



and the dual one is. 



(53) 



x*(0.,) = -^{e'^'--*-^-H^<5(0.,)}, (54) 

where v = e^"^ — 1. 

We thus give the conventional conjecture as follows, 
0,i, 

- {1 - (g- l)p}log{l -{q- l)p} -{q- l)p\ogp 
= 77 log g- (55) 



9-1 



,/3J{<5(0o+ioi)+'5(0o+'o2)+'5(0o+'o3)+'5(0o+'o4)} 



= q 



-4i) + z;2^,5(Zo^,^) 



(58) 



where i ^ j means the summation over different pairs 
among four bonds, and i ^ j ^ k expresses the summa- 
tion over all combinations of different three bonds among 
four bonds. In addition, the dual principal Boltzmann 
factor is given as. 



4 g-1 4 



00=0 !:=i 

„4 r . 4 



(59) 



The solutions are obtained as Pc — 0.079731 for q = i, 
Pc = 0.063097 for g = 4, and pc = 0.052467 for q = b. 



From these quantities, we rewrite Eq. (|23p as follows 



log 



{q + vf 



q + Av + v'^Y^6{Iqi,Ioj) + ^ 6{loi,loj,lok) + v'^S{loi,lo2,lo3,lo4:) 



= 21og(7. 



(60) 
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FIG. 11: Relationships between Y and ^^R' ^^'^ between Y* 
and Z^^-^\ The spin S' at the center of the star shape is of 
the summation for the star-triangle transformation. 



where the configurational average for the random vari- 
ables {loi} on the four bonds follows the distribution 



function 

.(1 



We obtain pi^^ = 0.0791462 for q = 3, 

Pc ' = 0.0626157 for <? = 4, and = 0.0520578 for 
q — b. These results are shown in Table HIl for compar- 
ison with those by the conventional conjecture and the 
existing result by a numerical estimation [28| . 



D. First Approximation for the Triangular Lattice 

We show the explicit calculation of the first approxi- 
mation by the improved technique for the ± J Ising model 
on the triangular lattice. We consider the cluster labeled 
by s = 1 with three up-pointing triangles as in Fig. [9l In 
this case, it is convenient to define the following quanti- 
ties, 



y{s.{J^,}) 



I J J 

p/3(,7oi + Jo2S-f J03S) 



(61) 



y*{s,{j,,}) 



S' = ±l i=l 

- j|n(2cosh/3Jo.)(e' 



+ S'Se 



,/3Jo3 



Se 



J|(2sinh/3Jo,)(e''^''' -S'e' 



(62) 



where the locations of the spins S' , S, and the sets of 
{Joi} are described in Fig. [TTJ The summation over 
S' corresponds to the star-triangle transformation. We 



calculate two partition functions for the cluster as, by use 
of two quantities Y and Y*, 

3 

= E I[nsoAJk,^J}) 

So=±l fe=l 



= 2e'^^'=i ^"-"^ cosh{/3^ (Jfe,o2 + Jfc.03)}, (63) 

k=l 

^Tr'(/5' {Jij}) 
3 



E 

S'o = ±l k=l 

( 3 



2^ n n cosh/3Jfc,o^ + n sinll/3Jfc,o^ 
3 / 2 

+ W ]^cosh/3Jfc_oj sinh/3Jfc^c 

k=l \i=l 
2 

+ ]^ sinh pjkfii cosh 13 Jkx 



;.03 



(64) 



Substituting these quantities into Eq. (|30p . we can 
obtain the first approximation by the improved tech- 
nique on the triangular lattice. The result is given as 

p^P = 0.835957, which is also listed in Table HI with the 

(21 

one Pc = 0.835985 by another approximation by the 
numerical manipulation of the s = 2 cluster in Fig. [9l 
Each estimation gives the precise locations of the mul- 

ticritical point on the hexagonal lattice p^P = 0.932611 

(2) 

and Pc = 0.932593, by the relation for the mutually 
dual pair [Toj . 

ir(pTR)+i/(PHEx) = 1, (65) 

where Hij)) is the binary entropy defined by 

H{p) = -_plog2p- (l-p)log2(l-_p). (66) 

Even on the triangular lattice, we can consider the 
Potts spin glass ^29j , and apply the improved technique. 
In this case, we can describe, by the improved technique, 
a more precise line consisting of the multicritical points 
on a two-dimensional plane of couplings representing two- 
body interactions and three-body interactions. 

VI. PERFORMANCE OF IMPROVEMENT 

A. Multicritical Point 

We here discuss the performance of the improved tech- 
nique. 

At first, we remark the predictions for the location 
of the multicritical point of the ±J Ising model on the 
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Type Conjecture Numerical result 



SQ ±J 




0.889972 [7, 8 


0.8905(5) [19] 




pi''- 


0.890725 


0.8906(2) [20, 21] 




(1,1) 

Pc 


= 0.890794 


0.8907(2) [22| 




pf' = 


0.890822 


0.8894(9) [23] 




p'y' 


= 0.890813 


0.8900(5) [24] 




7(0) _ 




0.89081(7) [25] 


SQ Gaussian 


1.021770 [7, 8] 1.02098(4) [21] 






1.021564 




TR ±J 




0.835806 [llj 


0.8355(5) [24] 




„(i) _ 

Pc — 


0.835956 






pf'^ 


0.835985 




HEX ±J 


Pf' = 


0.932704 [11] 


0.9325(5) [24] 




„(1) _ 
Pc — 


0.932611 








0.932593 




SQ Potts((j = 


3) = 


0.079731 \7, 8 


0.079-0.080 [28] 




pi'' = 


0.079146 




SQ Potts(g = 


4) p*"' = 


0.063097 [8] 






P« = 


0.062616 




SQ Potts(g = 


5) pi°' = 


0.052467 [8] 






P« = 


0.052058 





TABLE II: Several predictions for the location of the mul- 
ticritical point by the improved technique. SQ denotes the 
square lattice, TR expresses the triangular lattice, and HEX 
means the hexagonal lattice. 
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FIG. 12: (Color online) Phase diagram of the ± J Ising model 
on the square lattice by the s — 1 improved technique. The 
vertical axis is the temperature, and the horizontal axis is the 
probability for Jij — J > oi the ±J Ising model. The thick 
dashed line (blue) is by the conventional conjecture and the 
thick solid line (blue) is by the improved technique. The thin 
dashed line (red) is the Nishimori line. 



square lattice. All of the results by the improved tech- 
nique for the ± J Ising model on the square lattice indi- 
cate a higher value about Pc ~ 0.8908 than pc°'' ~ 0.8900 
by the conventional conjecture. As the size of the cluster 
increases, the prediction of Pc converges to some value 
about Pc ~ 0.8908. We need the precision to the fourth 
digit to conclude the conflict between pc ~ 0.8900p3l[2^ 
and Pc « 0.8908[il HO, [HI El El . For this purpose, 
the improved technique gives a satisfactory answer that 
the multicritical point is located at Pc ~ 0.8908. We can- 
not completely deny the possibility that the multicritical 
point locates at pc « 0.8900 as estimated in other stud- 
ies [H E3) because the improved technique does not 
give the exact solution. However the following discus- 
sions support our conclusion pc ~ 0.8908 from a different 
point of view. 



B. Phase Diagram 

The phase boundary can be predicted by the improved 
technique without the restriction of the Nishimori-line 
condition, similarly to the conventional conjecture. Un- 
fortunately the improved technique fails again to derive 
the precise phase boundary especially under the Nishi- 
mori line similarly to the case of the hierarchical lattice 
as in Fig. [T^l Nevertheless we find an improve- 



ment of estimations of the slope of the phase boundary 
at the critical point Tc of the non-random Ising model. 
We concentrate on the estimations of the slope and show 
the results below by use of not only (s = 1), (s = 2), 
{s = l,t — 1), and {s — 1, t = 2) clusters but also several 
ones as in Fig. 1131 The computing time of the order 

0(2^B ) is needed in general for the configurational av- 
erage over {Jij} in evaluation of the improved technique. 
However, the configurational average becomes much sim- 
pler, when we consider a calculation only around Tc to 
estimate the slope at Tc, at which at most a single bond 
becomes antiferromagnetic. Therefore we can deal with 
further approximations only to estimate the value of the 
slope at Tc by various clusters. The obtained values are 
listed in Table IIIII Two types of the approximations, 
by use of one clusters and by dividing the square lattice 
into two clusters, give different values but, in any cases, 
the increase of the size of the cluster shows convergence 
to the exact solution 3.20911 of the slope at Tc by Do- 
many [T|. Therefore it is considered that the improved 
technique gives a systematic way to derive the precise 
locations of the critical points in the region especially 
above the Nishimori line. 
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FIG. 13: Several approximations for estimations of the slope 
at Tc. The numbers of bonds in each pattern are 36 (s = 3), 
40 {t = 3), 64 (s = 4), and 60 {t = 4). 



C. Verticality and Reentrance 

We also examine the shape of the phase boundary at 
the multicritical point of the ± J Ising model. The follow- 
ing observation also supports the validity of the improved 
technique. Taking the derivative by (3 of Eq. (|40|) . we ob- 
tain the expression of the slope of the phase boundary at 



the multicritical point, here denoted by (3, 

1 dZ* 



d(3p 
dp 



E 



dp 



1 dZ 

¥^dp 



E 



1 dZ* I dZ 



(67) 

where we omit the approximation type s, and the ar- 
guments of Z and Z* for simplicity. One can find the 
quantities in the numerator in the right-hand side of this 
equation being equal to the exact internal energy on the 
Nishimori line 



1 dZ* 
^ ¥^~dfT 



y-^ 1 dZ 



= TVstanh/JpJ. (68) 



Hence the slope of the phase boundary at the multicrit- 
ical point should be vertical. In other words, the multi- 
critical point is located at a minimum value Pc on the 
phase boundary predicted by the improved technique. 
This statement is satisfied for any clusters. Therefore, 



Type Number Value of slope 





s = 




1 3.41421 




s = 1 




4 3.33658 




s = 2 




16 3.31272 




s = 3 




32 3.29352 




s = 4 




64 3.28161 




exact 




3.20911 fl5] 




Type 




Number Value of slope 




s = 




1 3.41421 


s 


= l,t = 


1 


4 + 12 3.31225 


s 


= l,i = 


2 


4 + 24 3.29414 


s 


= l,i = 


3 


4 + 40 3.28170 


a 


= l,i = 


4 


4 + 60 3.27287 




exact 




3.20911 fl5] 



TABLE III: Slope at the critical point Tc for the ±J Ising 
model on the square lattice. The top table gives the results 
by the clusters with many cross shapes denoted by a. The bot- 
tom table shows those by two clusters represented by (s,t). 
For comparison, we write the result by the conventional con- 
jecture denoted by s = 0. 



even if we consider the infinite size of the small lattice in 
which we can expect to obtain the exact answer, the ver- 
ticality at the multicritical point holds. The improved 
technique can give the consistent phase boundary with 
the predicted by the gauge symmetry 0, 0] . 

The second derivative yields a non-zero value of 
d^Pp/dP"^, which is proportional to the difference of the 
specific heat between ± J Ising models on two clusters, 



1 1 (dZV 
2^Z \dp) 




^ f dZ* 



V2^° dPp 



1 dZ 

W^dPp 



\ogZ 



(69) 



We do not have the exact value of the specific heat on 
the Nishimori line though we know an upper bound 
We cannot completely determine the shape of the phase 
boundary only by the improved technique. However we 
remark that the estimated values of the second deriva- 
tive become lower, if the cluster under consideration be- 
come larger as 0.956729 (s = 0), 0.753892 (s = 1), and 
0.737262 (s = 2). These positive values indicate that the 
phase boundary predicted by the improved technique be- 
come reentrant or vertical. The possibility of the phase 
boundary is indeed limited into whether vertical or reen- 
trant as rigorously shown by the gauge-symmetry argu- 
ment 
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D. Other Random Spin Systems 

If the improvement affects tlie predictions not only of 
the multicritical point but also for other critical points, 
we can apply the improved technique to random spin 
systems without gauge symmetry. The absence of the 
Nishimori line on the phase diagram does not permit us 
to rewrite Eq. (|23|) as Eq. (j40|) . which can give a relation 
between two entropies of the distribution of frustration. 
However the previous discussions by the duality for the 
cluster are applicable to various random spin models. We 
then give the critical points by the following equation 
even for random spin systems without the Nishimori line, 



N 



N^,'^ + 1 log 2, 



(70) 



where [F^ 



-)1 



is the free energy on the cluster, and 

Therefore 



represents that on the dual cluster 

our task to analytically derive the critical points in ran- 
dom spin systems is to estimate the difference between 
the two free energies on the cluster and its dual one. 

We apply the improved technique to the bond-diluted 
Ising model by the evaluation of Eq. ([70]) by use of the 
distribution function, 



P(Jy) = pS{J,j - J) + (1 - p)(5( J„). 



(71) 



The predicted phase boundary in Fig. [TJ] is not drasti- 
cally different from the one by the conventional conjec- 
ture [lOl ■ However we can find a significant difference by 
investigation of the values of the slope at similarly to 
the case of the ± J Ising model. We show the results for 
the slope at Tc for the bond-diluted Ising model in Table 
ITVl The estimated values for the square lattice shows 
convergence to the exact solution 1.32926 [lB|, similarly 
to the case for the ± J Ising model. 

In addition, we remark that the improved technique 
works very well for the critical points of the bond-diluted 
q-state Potts model, and g-state Villain model [sTl. From 
these points of view, we conclude that the improved tech- 
nique is also a systematic approach leading to the precise 
locations of the critical points in broader classes of the 
random spin systems. 




0.5 0.6 0.7 0.8 0.9 1.0 



FIG. 14: (Color online) Phase diagram of the bond-diluted 
Ising model on the square lattice by the s = 1 improved tech- 
nique. The vertical axis is the temperature, and the horizontal 
axis is the probability for Jij = J > 0. 



Type Number Value of slope 





s = 




1 1.34254 




s = 1 




4 1.33780 




s = 2 




16 1.33626 




s = 3 




32 1.33500 




s = 4 




64 1.33420 




exact 




1.32926 [16] 




Type 




Number Value of slope 




s = 




1 1.34254 


s 


= l,i = 


1 


4 + 12 1.33623 


s 


= l,i = 


2 


4 -f 24 1.33504 


s 


= l,t = 


3 


4-^40 1.33421 


s 


= l,i = 


4 


4 4-60 1.33362 




exact 




1.32926 [16] 



TABLE IV: Slope at Tc for the bond-diluted Ising model on 
the square lattice. For comparison, we write the result by the 
conventional conjecture denoted by s = 0. 



VII. CONCLUSION 

We proposed an improved technique applicable to the 
square, triangular, hexagonal lattices, and derived the 
precise locations of the multicritical points for the ± J 
Ising model, the Gaussian Ising model, and the Potts spin 
glass on the square lattice, as well as the ± J Ising model 
on the triangular lattice and the hexagonal lattice. This 
improved technique is still approximation for the location 



of the multicritical point. However we can enhance the 
precision of the approximation by the summation over 
spins in the cluster taken from the considered lattice, if 
we need the precise location of the critical points in a 
random spin system. This would open a way to analyti- 
cally derive the location of the critical points in random 
spin systems with very high precision. Unfortunately, in 
the low-temperature region under the Nishimori line, the 



16 



improved technique cannot give satisfactory answers yet. 
We solve this problem in the low-temperature region un- 
der the Nishimori line, and have to examine the validity 
of some hypotheses on the improved technique. 

In this paper, we restrict ourselves to the random spin 
systems in two-dimensional systems. However we can 
apply the duality to other dimensional systems. For ex- 
ample, the duality can transform the random-bond Ising 
model on the three-dimensional cubic lattice into the 
random-plaquette gauge model on the three-dimensional 
cubic lattice. The random-plaquette gauge model is an 
attractive one in terms of the quantum toric code (3^.[33|. 
An accuracy threshold to correct error of the quantum 
toric code corresponds to the location of the multicriti- 
cal point on the random-plaquette gauge model with the 
random couplings following the ± J distribution function 
on the three-dimensional cubic lattice. The conventional 
conjecture relates this threshold with the location of the 
multicritical point of the ±J Ising model on the three- 
dimensional cubic lattice Q. The improved technique 
also cannot directly derive such an accuracy threshold, 
but can make more precise relationship between the lo- 
cations of the multicritical points on the random-bond 
Ising model and the random-plaquette gauge model. 

As another direction of studies in the future, we should 



clarify the physical meaning of the equation consisting 
of the entropy of the distribution of frustration, which 
determines the location of the multicritical point. 
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